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We develop a nanoscale dynamical mean field theory (nano-DMFT) to deal with strong Coulomb 
interaction effects in physical systems that are intermediate in size between atoms and bulk mate- 
rials, taking into account the tunneling into nearby electrodes. Focusing on a simplified tree-like 
geometry, the usual DMFT loop simply stops when the finite lattice is fully covered, starting with an 
initial seed provided by the electronic environment at the boundary. To illustrate this nano-DMFT, 
we investigate the disappearance of the quasi-particle weight in a correlated nano-object near the 
Mott transition. In contrast to thermally-driven classical phase transitions, quantum effects lead to 
unexpected oscillations of the order parameter, related to the interference of coherent renormalized 
quasiparticles. This behavior also implies a spatially inhomogeneous Mott localization process at 
the nanoscale. 



The recent developments in the field of strongly corre- 
lated electronic systems are reaching today an interesting 
point of convergence. Coming from a macroscopic de- 
scription of bulk condensed matter, powerful many body 
techniques [H, 0] , possibly in combination with ab-initio 
methods, have been tailored to understand an ever grow- 
ing number of new compounds with surprising proper- 
ties Q . From a bottom-up perspective, the recent exper- 
imental progresses in designing artificial atoms (quan- 
tum dots) or manipulating single atoms and molecules 
that are probed by an electronic environment (electrodes) 
have been met with as well an impressive success [4| . 

How those two fields will merge together is a challeng- 
ing question raised by today's research, whether in semi- 
conducting devices with more complex lithographic pat- 
terns, or by the study of ultrasmall chunks of correlated 
materials in quantum transport experiments. Such de- 
vices ask however for the development of new theoretical 
methods as well, since the inclusion of correlation effects 
in transport calculations by modern numerical methods 
can only be achieved at the level of single atomic impu- 
rities, using Wilson's Numerical Renormalization Group 
method (NRG) 0, @, or simplified linear "molecules", 
by the Density Matrix Renormalization Group @, @] , but 
fails for higher-coordination entities. 

In the case of correlated materials in the thermody- 
namic limit, a philosophy based on a local effective de- 
scription has proved nevertheless very successful. Indeed 
the so-called dynamical mean field theory (DMFT) ex- 
actly maps (in the large coordination limit) a macro- 
scopic lattice problem of correlated fermions onto a sin- 
gle atomic level that feels the average effect of the other 
electrons moving nearby [H, In this case, the above- 
mentioned numerical techniques can be used to solve this 
simplified effective problem, and allow to capture quanti- 
tatively the effect of the Coulomb interaction on the low 
energy electronic properties: strong signatures indeed ap- 
pear in the electrical conductivity II, U3] as electrons 
start to localize near the Mott insulating state [lj]. 

The purpose of this paper is to develop a nano-DMFT 



approach for correlated nano-objects which are interme- 
diate in size between the atom and the solid. In the 
simplest case of tree-like structures, this can be easily 
formulated using the cavity construction and can be 
simply described by a standard self-consistency DMFT 
loop which is stopped at a finite number of iterations, 
related to the linear size of the nano object. This new 
insight will be applied both to the Ising and Hubbard 
models respectively, allowing to underline some impor- 
tant physical differences in finite-size many-body effects 
between the classical and the quantum world. As an 
illustration for the latter, we will examine how low en- 
ergy properties cross over from the single atomic impurity 
to the correlated solid as a function of the system size, 
focusing on the local quasi-particle weight, a quantity 
that characterizes the proportion of electrons that sur- 
vive the Coulomb repulsion at low energy, and can serve 
as an order parameter for the Mott transition. We find 
out a surprising oscillatory behavior of this quasiparticle 
residue as a function of physical parameters, which re- 
veals the interference of coherent quasiparticles in finite 
size systems. The Mott localization also happens to be 
non homogeneous at the mesoscale, with a proportion of 
sites considerably "hotter" than the others (in the sense 
of scattering properties). 

We begin by a general description of the cavity method 
for finite size lattices in the large coordination limit, see 
figure [TJ Although a derivation of the DMFT equations 
can be done in principle for an arbitrary geometry of the 
nano-object, we will consider for simplicity a Bethe lat- 
tice of linear size L (determined from its center) and co- 
ordination z\ this corresponds to the upper right panel of 
figure Q] for L = 2 and z — 4. Focusing on the properties 
at the center of the lattice (site 0), these can be com- 
puted from an effective local action involving the site 
only, provided one can obtain the dynamics at the nearby 
sites 1 in the lattice where the site has been removed 
(cavity), see the lower right panel of figure [T] This state- 
ment is exact in the limit where the coordination becomes 
large [l(| , but can be taken as a local mean field ansatz 
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FIG. 1: Correlated nano object of linear extension 2L + 1 = 5 
and coordination number 2 = 4, denned by a finite piece of 
either the square lattice (top left) or the Bethe lattice (top 
right). For the latter, which allows the simplest calculations, 
the nano-DMFT procedure is outlined by the two further 
stages of a cavity method given in the lower pictures. The 
external environment, to which the nano object couples, is 
pictured by the shaded area. The nano-DMFT becomes ex- 
act in the limit z — > oo keeping L finite. 
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FIG. 2: Average magnetization (So) at the center site for a 
small boundary field h/ J = 0.01 as a function of temperature 
T/J (abscissa) and number of sites L (ordinate). 



considering a Hubbard- Anderson model [8| for electrons 

H = -tJ2 d Ui„ - m E 4a. + u E 4aAAi 

ij,<7 icr i 

ka k(T 



otherwise. Similarly, the next shell can also be tackled 
using a action local to the sites 1, given by an effective 
field provided by the sites 2 (in the cavity where the sites 
of type and 1 have been removed) . This is repeated re- 
cursively until the boundary is reached, which amounts 
to an impurity coupled to a given external environment 
(lower left panel of figure [1]) . 

Let us put the reasoning backwards, focusing for 
the sake of clarity on the simple Ising model H = 
SiSj — hiSi, with a local external magnetic 
field that couples to the outer shell, = h ^ only 
(due to the tree structure of the lattice, all those sites 
are independent, so that an extra label is not required). 
These boundary spins are thus free, i.e. tfW = -hS L , 
which gives the thermal average (Sl) = tanh(/i/T). 
This corresponds to the lower left picture in figure [T] 
for L = 2. Now the physics of the sites L — 1 can 
be obtained by an effective impurity problem that reads 
H^^ 1 ) = —J(Sl)Sl-i (this is exact asz-^ oo), giving 
(Sl-i) = tanh( J(Sl)/T). This equation can be used 
recursively, until the center (site 0) is reached, provid- 
ing the results seen on figure O Besides the expected 
monotonous disappearance of the magnetization (So) a t 
large system size near the critical temperature T c / J = 1, 
we remark that the thermodynamic limit is only obtained 
for rather large system size (which increases with decreas- 
ing external magnetic field at the boundary). 

We now turn to the more interesting quantum case, 



with hopping t, chemical potential /i, local Coulomb re- 
pulsion U within the correlated lattice, while the hy- 
bridization V couples only the outer shell of sites L (an 
unnecessary extra index labelling the equivalent sites of 
type L has been omitted above) to an external electronic 
bath, which models the electrodes in a transport experi- 
ment. Following the previous philosophy, we start at the 
boundary (sites L) with single Anderson impurities (de- 
coupled from the other sites with i < L), so that the local 
Green function G^ L \u) is obtained by solving a local ef- 
fective action which is nothing but the single impurity 
Anderson model: 

S {L) = ~EH^n-M + AW(^ n )]4A CT 
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where f3 = 1/T, w„ = (2n + 1)ttT. Equation © pro- 
vides thus the boundary condition (and not an effective 
Weiss field at this first stage), similar to the magnetic 
field h applied at the boundary in the Ising model. Note 
that in the absence of an external electronic environment, 
A( L ) = 0, the boundary problem reduces to an atomic 
limit. The nano-DMFT for the full isolated molecule 
could then allow the calculation of its many-body level 
spectrum. 
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Now clearly the effective local action of the sites L — l 
takes a form similar to ([2]) , with the simple effective bath: 

^ L - 1 \uj n )=t 2 G^ L \iuj n ) (4) 

due to the Bethe lattice structure (this is obtained by 
integrating exactly the sites L in a cavity system con- 
stituted of sites L — l and L only, see the lower right 
panel of figure [TJ) . By repeating this procedure up to 
the center of the lattice, we see that the local Weiss field 
Green function at the center site in the Bethe lat- 
tice of linear size 2L + 1 is determined by the L th DMFT 
iteration loop with an initial condition (0) given by the 
external electronic environment. 

We now want to investigate some physical results that 
can be obtained from this nano-DMFT approach. For 
simplicity, we focus on the case where the probing elec- 
trodes are modelled by a constant density of states po, so 
that the external environment is characterized by a sin- 
gle width parameter T = irp V 2 , providing the starting 
Weiss field of the DMFT loop A iL ^(cj n ) = -iTsign(uj n ). 
This is then iterated up to a fixed number L of times, cor- 
responding to a given linear size L of the lattice, for each 
value of the Coulomb interaction U and temperature T 
(the chemical potential is chosen as p = U/2 to stay at 
particle-hole symmetry). In contrast, the usual DMFT 
procedure for a macroscopic system is done until a small 
and fixed convergence threshold is reached (including a 
weighting procedure that speeds up this convergence), 
and does not depend on the boundary condition (unless 
one considers the first-order metal/insulator transition 
regime where two solutions coexist, see below). 

Since the Hubbard model in its thermodynamics limit 
displays a quantum phase transition at increasing U from 
a correlated metal to a paramagnetic Mott insulator, we 
would like to look at a quantity which determines a sim- 
ple order parameter. Focusing on the site in the center 
of the small island, this is provided by the local quasi- 
particle residue Z^ at this site 0, defined as: 



l + idY,(°){iu)/du w = K ' 

introducing the local self-energy £(°)(iw) = ioj + fx - 
A(°)(w) - l/G^iiuj). For the infinite size Hubbard 
model, the quasiparticle weight is uniform, decreases 
from the value 1 at U = 0, and vanishes continuously 
at a critical strength U C 2/2t ~ 3. For U < U C 2, strongly 
renormalized quasiparticules with a reduced coherence 
time govern the low energy properties. In parallel to the 
stabilization of a finite magnetization in the Ising model 
with a non-zero boundary magnetic field, the coupling 
r of the outer sites to the electronic environment plays 
the role of a dispersion, and acts in stabilizing a con- 
ducting state, even in a finite size system. The quantity 
can thus be defined for all system sizes and interac- 
tions (even above U C 2 provided the thermodynamic limit 



is not reached). We are thus interested in the conver- 
gence of Z^ ' with the system size L by increasing the 
Coulomb interaction U . This is seen on figure [31 which is 
obtained by using the iterated perturbation theory (IPT) 
method 
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FIG. 3: Quasiparticule weight in the center of the 

molecule as a function of Coulomb interaction U/2t (abscissa) 
and system size L (ordinate), with external coupling to the 
bath T/2t — 0.001 and inverse temperature (3 = 800. 

Superficially, this plot is very similar to the magneti- 
zation curve of the finite size Ising model of figure [2] for 
any size L, the quasiparticle weight Z^> drops from 1 
as U increases; also similar to the Ising case, the conver- 
gence with L towards the thermodynamic result is slow 
when the coupling at the boundary to the environment 
is small. For this precise plot, this occurs roughly for 
L ~ 50, which for an actual Bethe lattice with connectiv- 
ity z = 4, amounts to a macroscopic system with N = 3 50 
sites, an unphysically too large size, due to the fact that 
the Bethe lattice is very hollow. In the more realistic 
case of a square lattice, one obtains a number of sites 
N = 100 2 , which is still very challenging for any direct 
numerical treatment of the complete many-body problem 
for such a molecule, providing a strong incentive in favor 
of the more economical nano-DMFT approach. 

There are however some important differences if one 
compares closely the classical result, Fig. [21 to the quan- 
tum one, Fig. [3] First, one sees that Z^> does not actu- 
ally continuously vanishes at L large, but presents a small 
jump at U c i/(2t) ~ 2.5. This is due to the existence in 
the range U c \ < U < U C 2 of both metallic and insu- 
lating solutions, to either of which the system converges 
depending on the initial condition. For the small cou- 
pling r considered here (corresponding to a tiny Kondo 
temperature for the initial outer sites L), one starts very 
close to an insulating seed, explaining the discontinuous 
character of this transition seen in figure [3J (the metallic 
branch is nevertheless correctly reached for U < U c \ and 
infinite DMFT iterations). One can check that the true 
vanishing of Z^ in the thermodynamic limit at U C 2 can 
be achieved by increasing the ratio Y /2t. 
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More interestingly, we discover by closer inspection of 
figure [3] some small oscillatory behavior of the quasi- 
particle weigth Z^ Q > as a function of U for intermedi- 
ate L sizes (beyond the transient regime of very small 
L). To display this effect more clearly, we plot on fig- 
ure 2] the relative weight to the thermodynamic limit 
SZ^ = Zf ] - These oscillations constitute a 

very surprising result that signals that localization (at 
site 0) is slightly non-monotonous with increasing inter- 
action, up to again very large system sizes. It can be 
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FIG. 4: Relative quasiparticle weight 8Z {0) = Zf^ - Z {0) 
for the same parameters as in figure [3] 

seen also that those oscillations are displaced when com- 
paring two adjacent sites, thus implying a non-uniform 
Mott localization process. Note that a similar behavior 
appears in the static prop erties of the Wigner localiza- 
tion in quantum dots [13(. The origin for these effects 
is the interplay between the interaction induced localiza- 
tion and the interferences within the nano object (due the 
presence of the boundary). This interpretation is consis- 
tent with the observation that these oscillations disap- 
pear both by increasing temperature, or by increasing 
the transmission at the boundary (i.e. making T com- 
parable to the hopping t). This illustrates nicely how 
correlation effects in intermediate size nano objects may 
lead to unexpected and complex physical behavior. 

We end the paper with some general remarks and 
desirable extensions of the nano-DMFT method. We 
would like first to emphasize that, despite some resem- 
blances, the nano-DMFT is not related to the cluster- 
DMFT scheme 14J. Indeed, in the latter one tries to 



compensate the sheer locality of DMFT by exactly solv- 
ing a finite size plaquette (with a total number N of sites 
certainly smaller than 10) embedded in a self-consistent 
medium that describes the effect of the remaining sites 
in an infinite lattice (it is a thus a cavity construction 
with a cluster of several sites taken as cavity object). 
The nano-DMFT allows however to reach very large pla- 
quette sizes, typical of big molecules, at the expense of 
neglecting some non-local aspects of the physics. Further 



extensions of the present work should focus on investigat- 
ing transport, as non homogeneities in the quasiparticle 
weight will also appear in the scattering rate for electron- 
electron interactions r~_ e cx (T/Z) 2 . This asks for more 
sophisticated solvers for the effective action © than the 
IPT method [H, e.g. with NRG [H, and the develop- 
ment of a consistent formalism for computing the tunnel- 
ing conductance. We stress finally that the hypothesis of 
a Bethe lattice is not central for the implementation of 
the nano-DMFT, and arbitrary "molecular" structures 
can be studied along the lines of Ref. [H, [HI, [l6| . This 
however needs to consider a fully iterated self-consistent 
nano-DMFT loop for the local Green functions, with as 
many local self-energies as the number of sites in the 
molecule. 

To conclude, we have extended the DMFT of macro- 
scopic correlated fermions models towards a nano-DMFT 
for generic nanoscale devices such as large molecules, 
quantum dot arrays, or nanostructures implanted with 
magnetic impurities. The method is naturally able to 
capture both strong electronic interactions and the cou- 
pling to an electronic environment. We have shown that 
these effects combine in mesoscopic systems to give in- 
teresting interplays in the Mott localization process that 
reveals the truly quantum nature of this phenomenon. 

The author thanks V. Dobrosavljevic, A. Lichtenstein 
and D. Ullmo for valuable discussions. 
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